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Abstract. Most of optical gravitational lenses recently discovered in the Sloan 
Digital Sky Survey Quasar Lens Search (SQLS) have two-images rather than four- 
images, in marked contrast to radio lenses for which the fraction of four-image lenses 
(quad fraction) is quite high. We revisit the quad fraction among optical lenses by 
taking the selection function of the SQLS into account. We find that the current 
observed quad fraction in the SQLS is indeed lower than, but consistent with, the 
prediction of our theoretical model. The low quad fraction among optical lenses, 
together with the high quad fraction among radio lenses, implies that the quasar optical 
luminosity function has a relatively shallow faint end slope. 
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1. Introduction 



Strongly lensed multiple quasars have been known to provide an unique probe of our 
universe. In particular, the point-source nature of quasars allows a simple statistical 
study from image multiplicities: Statistics of the number of multiple images provide 
constraints on the ellipticity and density profile of lens objects as well as the faint end 
luminosity function of source quasars [Il[2l[3j|IJ[5j[6l[71[El[9]. 

The statistics of image multiplicities have been done mainly using radio lenses. 
[5] adopted a radio lens sample of the Cosmic Lens All-Sky Survey (CLASS) JTUJ [11] 
to show that the fraction of four-image (quadruple) lenses is significantly higher than 
expected from a standard mass model of elliptical galaxies. [6] showed that the fraction 
of quadruple lenses in a statistical subsample of the CLASS is marginally consistent 
with what we expect from the observed galaxy population, but it still requires relatively 
large galaxy ellipticities. 

Recent large-scale optical surveys allow us to conduct complementary statistics 
using optical gravitational lenses. In particular, a large sample of quasars discovered in 
the Sloan Digital Sky Survey (SDSS) [12] is quite useful for a strong lens survey: Indeed, 
the SDSS Quasar Lens Search (SQLS) [13] has already discovered approximately 20 
new strongly lensed quasars (see, e.g., [H] and references therein), becoming the largest 
statistical sample of strongly lensed quasars. Interestingly, the fraction of four-image 
lenses (quad fraction) in the SQLS appears to be significantly lower than the CLASS. 
Only a few lenses among ~ 20 new SQLS lensed quasars are quadruple lenses, whereas 
nearly half of CLASS lenses were four (or more) image systems. 

In this paper, we revisit the quad fraction among optical gravitational lenses. We 
adjust the selection function to that of the SQLS and make a comprehensive prediction of 
the fraction of quadruple lenses. A particular emphasis is paid to whether the current low 
quad fraction in the SQLS is consistent with the observed galaxy properties. Throughout 
the paper we adopt A-dominated cosmology with the matter density VLm = 0.3 and the 
cosmological constant Qa = 0.7. 



2. Calculation 



2.1. Lensing Probabilities 



We assume that the mass distribution of galaxies can be approximated by an Singular 
Isothermal Ellipsiod (SIE). The scaled surface mass density of an SIE is given by 
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where e denotes the ellipticity. The Einstein radius 9e (for e = 0) is related with the 
galaxy velocity dispersion a by 
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with D\ s and D os being the angular diameter distance from lens to source and from 
observer to source, respectively. The normalization factor A(e) basically depends on 
the shape and viewing angle of galaxies: In this paper we assume that there are equal 
number of oblate and prolate galaxies and adopt the average of the two normalizations 
(see [6]). We find that with this normalization the Einstein radii are roughly equal for 
different ellipticities. 

It is expected that the quad fraction is mainly determined by the ellipticity. 
Although the external shear also produces the quadrupole moment in lens potentials, 
the effect is expected to be minor. For instance, the standard strength of external shear 
(median value of < 0.05) can cause notable changes in the quad fraction only for lens 
galaxies with e < 0.2 [5]. Therefore throughout the paper we neglect the external shear. 

We solve the lens equation using a public code lensmodel [15]. The lensing cross 

section a\ ens is computed by summing up source positions that yield multiple images with 

a weight of $(L//i) / yu/$(L) , where $(£) is the luminosity function of source quasars and 

/i is the magnification factor (see §2 .21 for which magnification factor we adopt). Lensing 

cross sections are derived for double and quadruple lenses separately. We compute the 

image separation for each event from the maximum separation between any image pairs. 

In computing the lensing probability, we impose a condition that the lensing galaxy 

should not be brighter than the source quasar, because the lens system may not be 

targeted for spectroscopy if the lensing galaxy dominates in the flux. We compute the 

galaxy luminosity from the velocity dispersion adopting an observed correlation [16J. 

Then the lensing probability of a source at z = z s becomes 

dpi r z * cdt 3 r da lcnSti dn ~ . 

- = J Q d Zl —(l + Zl ) J da^--6(e - 0)9(z gal - Vo ), (3) 

with dn/da being the velocity function of galaxies. The suffix i = 2 or 4 denote the 
number of images. 

In computing the lensing probability, we need to specify the lens galaxy population. 
Since strong lensing is mostly caused by early-type galaxies, particularly for strong 
lenses in the SQLS whose image separations are basically larger than 1", we only 
consider early-type galaxies. For the velocity function, we assume that of early-type 
galaxies derived from the SDSS (T71 [18]. More important for the quad fraction is the 
distribution of ellipticities. We adopt a Gaussian distribution with mean e = 0.3 and 
the dispersion a e = 0.16, which is consistent with observed ellipticity distributions of 
early-type galaxies [HI EUJ EH [22], [23], as a fiducial distribution. However we also vary 
the mean ellipticity, e, to see how the quad fraction depends on the ellipticity. 

2.2. Quasar Population and Selection Function 

The quasar luminosity function is another important element to make an accurate 
prediction of the quad fraction. As a fiducial luminosity function, we adopt that 
constrained from the combination of the SDSS and 2dF [24J: 
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Table 1. A current statistical sample of lensed quasars in the SQLS. iVi mg indicate 
the number of quasar images. 
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where a pure luminosity evolution with 

M*(z) = M*(0) - 2.5(k lZ + k 2 z 2 ) (5) 

is assumed. The parameters are (3^ = 3.31, (3\ = 1.45, = 1.83 x 10 _6 Mpc _3 mag _1 , 
M*(0) = —21.61, k\ = 1.39, and k 2 = —0.29. We convert rest-frame g-band magnitudes 
to observed i-band magnitudes using K-correction derived in |25j . 

The selection function of the SQLS was studied in detail in [13]. Since the statistical 
sample of lensed quasars is constructed from quasars with i < 19.1 and at 0.6 < z < 2.2, 
we restrict our calculation in this range. The magnification bias is computer assuming 
the image separation dependent magnification factor (see equation (14) in p3]). At 
9 > 1" the completeness is almost unity, but there is a small difference of completeness 
between double and quad lenses: To take this into account we include completeness 
4>i(6) in our calculation. In summary, we compute the numbers of double and quad 
lenses as 

n % = r dz s [ dM g <S>(M g )QD 2 os ^(l + z s ) 3 f" d6&(6)^, (6) 

where i\ im = 19.1 for the statistical lens sample of the SQLS. We have set the upper limit 
of the image separation to 3" since beyond the image separation the effect of surrounding 
dark matter becomes significant (see, e.g., |26j). The fraction of quadruple lenses is then 
computed as 

PQ = : • (7) 

n 2 + n 4 

2.3. Lensed Quasars in the SQLS 



The SQLS has already discovered about 20 new lensed quasars as well as several 
previously known lensed quasars. Although the statistical sample of lensed quasars 
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Figure 1. The fraction of quadruple lenses pq as a function of i-band limiting 
magnitude in m . Here we consider lensed quasars with redshifts 0.6 < z < 2.2, flux 
ratios > 10~ 05 , image separations 1" < 6 < 3", and lens galaxies fainter than 
the quasar components z ga i — i qso > 0. Dotted line indicate the limiting magnitude of 
SDSS quasars, i = 19.1. Left: From lower to upper solid lines, the faint end luminosity 
function of quasars [3\ is changed from 1.05 to 1.85. The mean ellipticity e is fixed to 
0.3. Right: From lower to upper solid lines, the mean ellipticity e is changed from 0.1 
to 0.5. The slope (3\ is fixed to 1.45. 



is still to be finalized, we use these lenses to make a tentative comparison with the 
theoretical expectation. To make a fair comparison with theory, we select a subsample 
of lenses by choosing lenses with redshifts 0.6 < z < 2.2, magnitudes i < 19.1, i- 
band flux ratios (for doubles) /, > 1CT 0,5 , image separations 1" < 9 < 3", and lens 
galaxies fainter than the quasar components 2 ga i — i qso > 0. Currently we have 13 lensed 
quasars that meet these conditions, which are summarized in Table [TJ Among these 13 
lenses only two are quadruple lenses, thus the observed quad fraction for the flux limit 
i lim = 19.1 is p Q = 2/13 ~ 0.154. 



Before comparing our calculation with the observed quad fraction, we see how it 
depends on parameters. Among others, the most important parameter is the ellipticity. 
Another important element that determines the quad fraction is the shape of the quasar 
luminosity function. In particular the faint end slope f3\ still contains large errors because 
current large-scale surveys are not deep enough to fully explore the faint end luminosity 
function. For instance, [37] and [38] adopted the 2dF quasar sample to derive the faint 
end slopes of f3\ = 1.58 and 1.09, respectively A survey of faint quasars conducted by [39] 




3. Result 
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Figure 2. The quad fraction in our fiducial model (e = 0.3, f3\ = 1.45; shown by 
a solid line) is compared with observed fractions in the SQLS (filled triangles with 
errorbars). The errors indicate 68% error estimated assuming the Poisson distribution 
for the numbers of double and quad lenses. See table Q] for the lens sample we use. 
Note that the data points are not independent but rather correlated in the sense that 
lenses used to plot at each i\ im are included in computing data points at larg er &iim as 
well. 

suggests that the faint end slope could be (3\ = 1.25, shallower than our fiducial value. 
Other uncertainties, such as cosmological parameters, the velocity function of galaxies, 
and the number of source quasars, affect the number of double and quad lenses roughly 
similarly, thus they hardly change the fraction of quad lenses. We find that the effect 
of changing the prolate/oblate fraction is not large, affecting the quad fraction only by 
a few percent. Therefore, in figure [1] we plot the quad fraction as a function of the 
limiting magnitude in m changing these two important parameters, the mean ellipticity e 
and the faint end slope (3\. First, the quad fraction decreases as the limiting magnitude 
increases. Larger magnifications of quads than doubles indicate that the quad fraction 
is a strong function of magnification bias such that larger magnification bias results in 
larger quad fraction, which explain the decrease of the quad fraction with increasing 
iw m . As expected, the quad fraction is quite sensitive to the ellipticity and the faint end 
slope of the quasar luminosity function. 

Next we compare the quad fraction in our fiducial model with the observed fraction 
in the SQLS. Figure [2] shows both the theoretical and observed quad fractions as a 
function of the limiting magnitude. We find that the observed quad fraction is indeed 
lower than the theoretical prediction. For instance, at i\\ m = 19.1 the quad fraction in 
our model is pq = 0.273 that is larger than the observation, pq = 0.154. However, by 
taking the large errorbar of the observed fraction due to the small number of lenses, we 
conclude that the observed quad fraction is consistent with the theoretical expectation. 
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Figure 3. Probability of our model producing quad lenses equal or fewer than iV qua d 
in a sample of 13 lenses, computed from our model prediction of the quad fraction for 
'him = 19.1, Pq = 0.273. In addition to our fiducial model we also plot the probability 
for e = 0.4 (pq — 0.408) that better reproduces the high quad fraction in the CLASS. 
The observed number of quads in the SQLS, iVq Ua d = 2, is indicated by a vertical 
dotted line. 



In figure [3] we plot the probability that our theoretical model produces < iV qua H 
quad lenses in a sample of 13 lenses. Note that in observation there are iV quad = 2 quad 
lenses (see tabled]). In our fiducial model the probability is ~ 0.27, which is low but 
acceptable. On the other hand, if we increase the mean ellipticity to e = 0.4, which 
is roughly the best-fit value for the observed quad fraction in the CLASS (see j6]), the 
probability reduces to ~ 0.05. Therefore with such large-ellipticity model it is difficult 
to account for the low quad fraction observed in the SQLS. 

Finally we check the dependence of the likelihood for i\q Ua d < 2 in a sample of 
13 lenses on the faint end slope j3\ and the mean ellipticity e in figure HI As expected 
from figure [U, the probability depends sensitively on these parameters. For instance, 
by decreasing the faint end slope to j3\ = 1.25, which is preferred by a spectroscopic 
survey of faint quasars [39], the probability is increased to ~ 0.35. Changing the mean 
ellipticity to 0.2 enhances the probability to ~ 0.61, making the observed low quad 
fraction quite reasonable. 

4. Summary and Discussion 

In this paper, we have studies the fraction of four-image lenses among optical 
gravitational lenses. We have paid a particular emphasis to whether the low quad 
fraction observed in the SQLS is consistent with the standard theoretical prediction. 
In order to make a fair comparison, we have taken account of the selection function 
and source population in predicting the quad fraction. We find that the observed quad 
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Figure 4. Probability of our model producing quad lenses equal or fewer than the 
observed case, 2V qua d = 2, is plotted as a function of f3i (left) or e (right). The fiducial 
values are shown by vertical dotted lines. 

fraction in the SQLS, pq = 2/13 ~ 0.154, is indeed lower than the prediction of our 
fiducial model, pq = 0.273, but is consistent given the large Poisson error of the observed 
quad fraction. 

We can lower the expected quad fraction by either making the faint end slope 
of the quasar luminosity function shallower or decreasing the mean ellipticity of lens 
galaxies. However, lowering the ellipticity decreases both optical and radio quad 
fractions, therefore such models have difficulty in explaining the high quad fraction 
among CLASS lenses. For instance, from the CLASS lens sample [6] derived 68% lower 
limit of the mean ellipticity to 0.28 which is marginally consistent with our fiducial 
model, e = 0.3. Therefore, one way to explain both the high quad fraction among radio 
lenses and low quad fraction among optical lenses is to consider a shallow faint end slope 
of the quasar optical luminosity function while keeping the mean ellipticity relatively 
high. 

A caveat is that the SQLS is still ongoing and the lens sample is not yet finalized. 
We should use a final, larger lens sample of the SQLS to draw a more robust conclusion 
from the quad fraction. The final statistical sample is expected to contain roughly twice 
the number of lenses we used in this paper, thus the statistical error should be reduced 
significantly. 

Acknowledgments 

I thank Naohisa Inada and Kyu-Hyun Chae for discussions, and an anonymous referee 
for many suggestions. This work was supported in part by the Department of Energy 



Is there a quad problem among optical gravitational lenses? 



9 



contract DE-AC02-76SF00515. 



References 

Blandford R D and Kochanek C S 1987 Astrophys. J. 321 658 
Kormann R, Schneider P and Bartelmann M 1994 Astron. Astroph. 284 285 
Kochanek C S 1996 Astrophys. J. 473 595 

Keeton C R, Kochanek C S and Seljak U 1997 Astrophys. J. 482 604 
Rusin D and Tegmark M 2001 Astrophys. J. 553 709 
Chae K H 2003 Mon. Not. Roy. Astro. Soc. 346 746 
Cohn J D and Kochanek C S 2004 Astrophys. J. 608 25 
Oguri M and Keeton C R 2004 Astrophys. J. 610 663 
Huterer D, Keeton C R and Ma C P 2005 Astrophys. J. 624 34 
Myers S T et al. 2003 Mon. Not. Roy. Astro. Soc. 341 1 
Browne I W A et al. 2003 Mon. Not. Roy. Astro. Soc. 341 13 
York D G et al. 2000 Astron. J. 120 1579 
Oguri M et al. 2006 Astron. J. 132 999 
Inada N et al. 2007 Astron. J. 133 206 
Keeton C R 2001 |astro-ph/0102340| 
Bernardi M et al. 2003 Astron. J. 125 1849 
Choi Y Y, Park C and Vogeley M S 2007 Astrophys. J. 658 884 
Chae K H 2007 Astrophys. J. 658 L71 

Bender R, Surma P, Doebereiner S, Moellenhoff C and Madejsky R 1989 Astron. Astrophys. 217 
35 

Saglia R P, Bender R and Dressier A 1993 Astron. Astrophys. 279 75 
Jorgensen I, Franx M and Kjaergaard P 1995 Mon. Not. Roy. Astro. Soc. 273 1097 
Rest A, van den Bosch F C, Jaffe W, Tran H, Tsvetanov Z, Ford H C, Davies J and Schafer J 

2001 Astron. J. 121 2431 
Sheth R K et al. 2003 Astrophys. J. 594 225 
Richards G T et al. 2005 Mon. Not. Roy. Astro. Soc. 360 839 
Richards G T et al. 2006 Astron. J. 131 2766 
Oguri M 2006 Mon. Not. Roy. Astro. Soc. 367 1241 
Inada N et al. 2005 Astron. J. 130 1967 
Inada N et al. 2006 Astron. J. 131 1934 

Oscoz A, Serra-Ricart M, Mediavilla E, Buitrago J and Goicoechea L J 1997 Astrophys. J. 491 L7 
Inada N et al. 2003a Astron. J. 126 666 

Schechter P L, Gregg M D, Becker R H, Helfand D J and White R L 1998 Astron. J. 115 1371 
Oguri M et al. 2005 Astrophys. J. 622 106 
Pindor B et al. 2006 Astron. J. 131 41 
Weymann R J et al. 1980 Nature 285 641 
Inada N et al. 2003b submitted 

Oguri M et al. 2004 Publ. Astron. Soc. Japan 56 399 

Boyle B J, Shanks T, Croom S M, Smith R J, Miller L, Loaring N and Heymans C 2000 Mon. 

Not. Roy. Astro. Soc. 317 1014 
Croom S M, Smith R J, Boyle B J, Shanks T, Miller L, Outram P J and Loaring N S 2004 Mon. 

Not. Roy. Astro. Soc. 349 1397 
Jiang L et al. 2006 Astron. J. 131 2788 



